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ABSTRACT 



Context. Swift detectors are found to be more sensitive to long-soft bursts than pre-Swift missions. This may largely bias the distri- 
bution of durations and then classification of gamma-ray bursts. 

Aims. In this paper, we systematically investigate the duration distribution of gamma-ray bursts in the Swift era via comparison with 
that of pre-Swift bursts. 

Methods. For the purpose of this study, statistical methods such as the K-S test and linear/non-linear fitting analysis have been used 
to examine the duration properties of Swift bursts in both observer and source frames. 

Results. For 95 GRBs with known redshift, we show that two-lognormal distributions of duration are clearly divided at Tgo -2 s. 
The intrinsic durations also show a bimodal distribution but shift systematically toward the smaller value and the distribution exhibits 
a narrower width compared with the observed one. Swift long bursts exhibit a wider duration dynamic range in both observer and 
source frame in comparison to pre-Swift long bursts. 

Conclusions. We find that Swift bursts and pre-Swift ones can share the same criterion of classification in terms of duration at 2 
seconds, although both monitors have large difference with respect to sensitivity of a given energy band. 
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1. Introduction 

Cosmic gamma-ray bursts (GRBs) are the most violent explo- 
sions occurring at cosmological distances in the universe. When 
a GRJB takes place, satellites can monitor its temporal variability 
in the y-ray energy band. The duration of the burst, 7*90, is de- 
fined as the time interval in which the integrated photon counts 
increase from 5% to 95% of the total counts. Based on an anal- 
ysis of durations using initial BASTE data, Kouveliotou et al. 
(1993) divided GRBs into two classes, i.e., long GRBs (LGRBs) 
with T 90 > 2 s and short GRBs (SGRBs) with T 90 < 2 s. The 
dichotomy has been justified by subsequent investigations (e.g., 
Mao, Narayan & Piran 1994; Katz & Canel 1996; Meegan et 
al.1996; Paciesas et al. 1999; Fishman 1999). In fact, the best 
parameters of a two-component lognormal fit to the distribution 
data were first obtained by McBreen et al. (1994). This fit has 
been supported by the current BATSE data with peak flux infor- 
mation (Horvath 2002; Nakar 2007). 

So far, much evidence showing the difference between two 
classes has been discovered and presented (see Zhang 2006 for 
review). For example, the spectra of LGRBs are softer than that 
of SGRBs in general. Besides, the pulse profiles of SGRBs are 
on average more symmetric than those of LGRBs (Zhang & Xie 
2007). The current Swift observations show that LGRBs have 
their median cosmological redshift z m ~ 2.0 higher than that 
of SGRBs z m ~ 0.4. All these differences suggest that both 
LGRBs and SGRBs, most likely, are distinct physical phenom- 
ena and produced due to model-independent emission engines 
(e.g., Balazs et al. 2003). However, it is not clear what causes 



these differences, especially the origin of bimodal Tgg distribu- 
tion. 

Koshut et al. (1996) pointed out that the observed duration 
distribution may vary with instruments. It is therefore neces- 
sary to investigate if there exists a new GRB class and/or what 
physical factors produce such properties (Gehrels et al. 2004). 
Naturally, we focus our study on the related issues of GRB clas- 
sification using the updated Swift data^^B 



2. Observation and data selection 

The higher sensitivity and angular resolution of Swift make it su- 
perior to previous space telescopes (e.g. BATSE, BeppoSAX and 
HETE-2), accurate follow-up observations have further added to 
its advantage (Gehrels et al. 2004). It can detect on an average 
about 2 GRBs per week within a 2 sr field of view, which is 
about two times more than that of pre-Swift missions (Meszaios 
2006). However, the detection rate of SGRBs to total GRBs is 
much lower (~ 8%) than the rate by BATSE (~ 18%), which 
is attributed to both their different energy responses and the rel- 
atively high spectral hardness of SGRBs (e.g., Band 2006a, b; 
Gehrels et al. 2007). On the other hand, the lower sensitivity to 
short duration bursts of Swift, relative to BATSE, makes it accu- 
mulate relatively lower counts, comparable with the number of 
background counts (Band 2006a). The effect of the instrument 
may cause the detection rate of SGRBs to be somewhat under- 
estimated. 
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Table 1. GRB samples of known durations (Tgo), redshifts (z) 
and/or peak energies (E p ) 



Sample 


No. 


Class 


Parameters 




Ref. 


sl + 


75 


LGRBs 


T90; z 




l; l; - 


s2- f 


20 


SGRBs 


T90', z 




l; l; - 


s3+ 


44 


LGRBs 


Tg ; z; E p 




1; 1;4 


s4' f 


11 


SGRBs 


T 90 ; z; E p 




1; 1;4 


s5* 


48 


LGRBs 


Tg ; z 


2, 


3; 5;- 


s6» 


18 


LGRBs 


Tg ; z 


2, 


3; 5;- 



Note: The references are given in order for Tg , z and E p . Data sets: f 
Swift; * pre-Swift; » BATSE. Of the 20 SGRBs, GRB 050709, is 
adopted from HETE-2. 

Reference: 1. Swift public webpage; 2. Ghirlanda, Ghisellini & Lazzati 
2004; 3. GRBlog website; 4. Butler et al. 2007; 5. Friedman & Bloom 
2005. 



In order to study the intrinsic properties of GRBs, we se- 
lected six data sets, namely si - s6, as listed in Table 1. As of 
2007 July 1, Swift has detected 75 LGRBs (si) and 20 SGRBs 
(s2) with known duration and redshift. The E p in the observed 
vF v energy spectra are chosen to characterize the spectral hard- 
ness relations with duration. Here, 44 GRBs from si have also 
the available E p values and constitute our sample s3. Out of the 
s2, only 1 1 sources have the measured E p and are employed to 
build the sample s4. Unfortunately, the remaining 9 bursts in s2 
(i.e., GRB 050202, GRB 050906, GRB 050925, GRB 051 105 A, 
GRB 051210, GRB 060313, GRB 070209, GRB 070406 and 
GRB 070429B) do not have the measured redshifts. For these 
we assigned a redshift value of z=0.5 to the 9 bursts, approaching 
the median redshift of z=0.4, as assumed by Norris & Bonnell 
(2006). These sources were included in the present study to im- 
prove the statistics. However, we found that the choice between 
z=0.5 and z=0.4 contributes only a small relative error of ~ 0.07, 
implying that the final results are not sensitive to the above as- 
sumed redshift values. In our fifth sample set, s5, we include 
48 pre-Swift LGRBs whose z and Tqq are already measured, in 
which 18 sources, less than half (~ 38%) of the 48 pre-Swift 
bursts, are detected by the BATSE mission and constitute our 
sample s6. 

3. results 

3. 1. Distributions of durations in the Swift era 

To check whether the distributions of Tgo are significantly differ- 
ent, indicative of a dependence on the on-board instruments, we 
first obtained the distributions of Swift bursts in both observer 
and source frames. Then, we compare the durations of LGRBs 
between Swift and pre-Swift/BATSE missions. 

3.1.1. Observed Tg Q distribution 

The accuracy of Tgo measurements is in principle affected by 
several factors including the identification of the time interval of 
a burst, the sensitivity of instrument, background modeling, the 
time resolution of the data, and the detailed shape of the burst 
time profile, etc. Figure 1 shows the Tgo distribution for the 95 
Swift GRBs, which include si and s2 samples (see Table 1). The 
best fit with a two-lognormal function gives the center values 
<T 9 oj,i = s and T 90 , p2 = 42.83^* s ) and the widths 

(wi = 19.05^4-60 s and w 2 = 18.20+^ s) with the reduced Chi- 




-2-10123 



iogT m „ (s) 

3 90, obs v ' 

Fig. 1. Bimodal distribution of durations for the 95 GRBs (si 
and s2; histogram) and two-component lognormal fit to the data 
(solid line). The GRBs are divided into two classes at Tgo ^ 2 s 
(vertical line): LGRBs (dotted line) and SGRBs (dashed line). 



square ^- 2 /dof = 0.67, which are roughly consistent with those 
calculated from the BATSE data (McBreen et al. 1994; Meegan 
et al. 1996; Paciesas et al. 1999; Horvath 2002; Nakar 2007). 
Note that the number of objects classified as belonging to the 
two lognormals to itself has been allowed to be a free parameter. 
The superposed function has a minimum around 2 s as found by 
Kouveliotou et al. (1993), indicating that the Swift sources are 
also divided into two classes, SGRBs and LGRBs, although the 
Swift is more sensitive to long soft bursts than the BATSE (Band 
2006a,b; Gehrels et al. 2007). This is an interesting result since 
the different detectors with diverse bandwidth sensitivity do not 
indeed affect the classification of GRBs in terms of duration. 

3.1.2. Intrinsic r 90 distribution 

One of the great Swift progresses is the increase in number of 
higher redshift sources. It is known that the median redshift of 
Swift LGRBs, z m ~ 2, is roughly two times larger than that 
of pre-Swift. As pointed out by previous authors (Bagoly et al. 
2006; Jakobsson et al. 2006), the difference of redshift distri- 
butions between the two samples is statistically significant. This 
may lead to an evident discrepancy between the two intrinsic Tgo 
distributions. We therefore utilize the Swift sources to explore 
this possibility. 

The potential spectrum evolution, as a result of cosmological 
redshift, can usually cause high energy y photons to shift into or 
out of the finitely sensitive bandwidth of a given detector. Note 
that here we had neglected the effect of a burst's spectrum itself 
softening with time (Ford et al. 1995). In this case, the trans- 
formation of Tyo from observer frame to source frame is gener- 
ally expressed as Tgoj nt =Tgol{\ + z) w , in which lo = 0.6 or 1, 
corresponding to energy stretching or not (Fenimore & Bloom 
1995; Meszaros & Meszaros 1995; Meszaros & Meszaros 1996; 
Horvath, Meszaros & Meszaros 1996; Balazs et al. 2003). Here, 
we consider a simple case of u> — 1 throughout this work, which 
the intrinsic duration is given by Tgoj nt =Tgo/(l + z). We calcu- 
late the intrinsic duration (T^int) distribution for the 95 GRBs 
and compare it with the observed one. As shown in Figure 2, the 
Tgoj nt has a bimodal distribution and is significantly shifted to- 
ward shorter durations than the observed one. The best fit with a 
two-lognormal function gives two centers (Tgo, p \ = 0.13*91? s, 
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Fig. 2. Comparison between the observed (solid line the same as 
Fig. 1) and the intrinsic (dashed-dotted line) Tgrj distributions. 



r 9 o, p2 = s ) an d two widths (wi = 10.96+^jf s and 

vv 2 = 17.38^ g3 s) with^ 2 /dof - 0.92, indicating that the distri- 
bution of TgQ,mt i s indeed bimodal but systematically narrower 
and shifted towards low value of durations in comparison to 
the observed one. A Kolmogorov-Smirnov (K-S) test returns the 
statistic D = 0.27 with a probability of P — 0.001, suggest- 
ing the intrinsic and observed duration distributions are drawn 
from the different parent populations. This result is well consis- 
tent with the theoretical prediction by Bromm & Loeb (2002), 
where they assumed the formation of all GRBs tightly follows 
the cosmic star formation history. However, some instrument- 
dependent factors such as the selection effects and discrepancy 
between different detectors can distort the resultant parent dis- 
tribution of Ti)oj nl in a sense. Once the disturbances are really 
reduced from the observed distribution, the remainder is the true 
parent distribution correlated with some physical predictions. In 
this case, our current result can offer some corroborative statisti- 
cal evidence that the rate of GRBs may really trace the star for- 
mation history, partly because the redshifts of SGRBs had not 
been measured before Swift. 

3.1 .3. Tg Q of LGRBs: Swift vs. pre-Swift 

Tgo may be wrongly estimated because of the factors pointed out 
in section 3.1.1. In particular, the measured duration has been 
found to be instrument-dependent (e.g., Koshut et al. 1996). We 
examined the dependence of Tc^ on the instruments by contrast- 
ing the normalized distributions of different detectors in Figure 
3. For this comparison, three samples including si, s5, and s6 
are applied. 

A lognormal fit to the six distributions gives the best pa- 
rameters listed in Table 2. We find the Swift sample obviously 
exhibits a wider duration distribution than the pre-Swift and 
BATSE did in both observer and source frames, while the widths 
of pre-Swift and BATSE distributions of duration are robustly 
equal. This implies that Swift can detect GRBs in a wider dy- 
namic range of Tgo, i.e., the fraction of longer and shorter Swift 
LGRBs are significantly higher than that of pre-Swift/BATSE 
LGRBs, except the increased detection rate of LGRBs due to 
the higher sensitivity of BAT loaded on Swift to long soft bursts 
(Band 2006 a,b). Considering the fact that Swift bursts can still 
be separated at Tgo = 2 s, we predict that the width of Tgo dis- 
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Fig. 3. The normalized duration distributions for the Swift (si; 
dashed line), pre-Swift (s5; solid line) and BATSE (s6; dotted 
line) LGRBs in the observer (upper panel) and the source (lower 
panel) frames. 



tribution of Swift SGRBs could be equally wider in comparison 
to pre-Swift detectors. Furthermore, we see from Table 2 that 
the distribution centers of Tgo,„b s are very close to each other, 
while for the Tgoj„t distribution, the centers of the pre-Swift and 
BATSE bursts are much more close but significantly larger than 
that of Swift bursts. It is interesting to note that the K-S test to 
the pre-Swift and BATSE bursts provides a very large probabil- 
ity of P ~ 0.99 in both observer and source frames showing that 
the s5 and s6 samples are consistent with being drawn from the 
same parent distribution. As for the Swift Tgoj nt distribution, the 
obvious decrease of the center value is attributable to the relative 
increase of the fraction of high redshift sources, which could be 
caused by the ubiquitous threshold effect of different instruments 
(e.g., Bromm & Loeb 2002; Band 2006a,b). 

3.2. Spectral hardness relation 

In theory, the peak energy E p can be used to reflect the photon 
components, similar to hardness ratio between different energy 
channels. The intrinsic peak energy (E pi ) in the source frame is 
related with the E p in the observer frame by E p j - E p (l + z). To 
verify the relation of spectral hardness with duration, we made 
plots of peak energy versus Tgo t0 b s for the samples s3 and s4 as 
shown in Figure Ab and Ac. A linear correlation analysis gives 
the coefficient of r — -0.27 with probability P - 0.07 for Figure 
Ab and r = 0.004 with P = 0.98 for Figure Ac, indicating that the 
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Table 2. The best-fit parameters to the three data sets 







logrgo.ois 






logT 90 j„, 




Sample 




W 






w 




si 


1.57±0.04 


1.31+0.08 


2.4 


1.06+0.04 


1.23+0.09 


2.8 


s5 


1.63±0.04 


0.90±0.09 


2.5 


1.32±0.03 


0.89+0.05 


0.9 


s6 


1.62±0.04 


0.76±0.07 


0.5 


1.29±0.03 


0.92±0.07 


0.3 



Note: The fitted parameters for the distribution center (jj) and the width (vv) are given in a logarithmic scale. 




Fig. 4. Relations of Tgo with E p (panel b) and E p j (panel c). Panels a and d show the distributions of the E p and E p j in logarithmic 
timescale, respectively. 



tendency of SGRBs relatively harder than LGRBs exists in the 
observer frame and almost disappears in the source frame. The 
former accords with the previous findings, for example, using 
the BATSE data (e.g., Kouveliotou et al. 1993) or the Swift data 
(e.g., Band 2006a). 

We also compare the peak energy of LGRBs and SGRBs in 
both observer (Figure 4a) and source frame (Figure 4d). It needs 
to point out that more cosmological redshift correction is re- 
quired to give LGRBs when we convert the physical variables 
from the observer frame to the source frame. As a result, the 
trend of the relative spectral softness for LGRBs (or the relative 
hardness for SGRBs) weakens more or less. For our samples, 
the long and short bursts have the median values of 74^H° keV 
and 398^20 keV, respectively, in the observer frame, and 302*$m 
keV and 617^^ keV, respectively, in the source frame. A K-S 
test returns the different probabilities of P — 0.01 for Figure 4a 
and P = 0.15 for Figure Ad, which hints that the E p j of long and 
short bursts may be drawn from the same parent population. 

4. Discussion 

The duration distribution from the third BATSE catalog is sug- 
gested to have a three-component lognormal form (Horvath 
1998; Mukherjee 1998; Hakkila et al. 2000). Similarly, the du- 
ration distribution of the current 4B catalogs exhibits the possi- 
ble existence of an intermediate group. Further studies show that 
the third group is either the excess of SGRBs with low fluence 
(Hakkila et al. 2003) or the softest LGRBs (Horvath et al. 2006). 
It is still controversial whether there exists the third class or not, 
due to the lack of a physical explanation (Horvath et al. 2006; 
see however, Chattopadhyay et al. 2007). The third group might 
be caused by an instrumental bias which reduces the durations 

4 http://www.batse.msfc.nasa.gov/batse/grb/duration/ 



of faint LGRBs (e.g., Hakkila et al. 2000; Hakkila et al. 2003), 
suggesting their presence is not physical but phenomenological. 
Therefore, the bimodality of the duration distribution is widely 
accepted by most people today. 

Basically, what causes the bimodality is an interesting but 
unsolved problem owing to the absence of direct observational 
evidence, although several works have tried to explain the bi- 
modality within different scenarios. It was suggested that the 
different spin axes of millisecond pulsars can interpret the two 
GRB classes (Usov 1992; Yi & Blackman 1998). Subsequently, 
Huang et al. (2003) studied the neutron star kick model (Dar & 
Plaga 1999) in detail and proved that this model can successfully 
account for the two-lognormal distributions if the central engine 
has a neutron star of high kick velocity larger than ~ 1000 km 
s -1 . Yamazaki, Ioka & Nakamura (2004) put forward the so- 
called unified model consisting of multiple subjets within an 
inhomogeneous main jet. Using this model, Toma, Yamazaki 
& Nakamura (2005) explained that the bimodality originates 
from discrete emitters in the main jet. They also predicted that 
the two kinds of bursts should have the same origin, i.e., su- 
pernovae, instead of the leading models, which predict that the 
LGRBs and SGRBs are produced due to the core collapse of a 
massive star and the merger of double compact objects, respec- 
tively (see e.g., Cheng & Lu 2001; Zhang & Meszaros 2004; 
Piran 2005; Meszaros 2006 and Lee & Ramirez-Ruiz 2007 for 
reviews). However, further investigations on the unified model 
showed that the bimodal distribution could be reproduced only 
for some special parameters (Janiuk et al. 2006). 

Recently, from an independent analysis of distinct 
timescales, Zhang et al. (2007) suggested that LGRBs oc- 
cur at larger distances from the central engine while SGRBs 
at smaller distances, i.e., two distinct y-ray emission regions 
may result in two different properties of GRBs, including the 
varieties of pulse profiles, as mentioned above. In this study, 
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they also pointed out the fact that LGRBs usually have long 
positive spectral lag (e.g., Norris, Marani & Bonnell 2000; 
Daigne & Mochkovitch 2003; Chen et al. 2005) and SGRBs 
have negligible lag (Norris & Bonnell 2006; Zhang et al. 
2006) can be naturally explained under the assumption that 
the curvature effect is a main contributor to the spectral lags 
(Ryde 2005; Shen, Song & Li 2005). The predication of distinct 
emission regions for different GRB classes and/or the nature 
of the bimodality still need to be clarified with more accurate 
observations in the future. 



5. Conclusions 

- We find that the observed durations of Swift GRBs have two- 
lognormal distributions divided clearly at Tw 2 s. This im- 
plies that the classification in terms of duration is unchanged 
from pre-Swift to Swift era. 

- The intrinsic durations also show a bimodal distribution but 
shifted systematically toward the smaller value and the dis- 
tribution exhibits a narrower width relative to the observed 
distribution. 

- The Tgo,int distributions of Swift and pre-Swift/BATSE 
LGRBs are significantly different in the width and center val- 
ues. 

- Swift bursts have a wider duration dynamic range than pre- 
Swift and BATSE bursts. 

- The spectra of the SGRBs are predominantly harder than the 
LGRBs, confirming the previous results from the BATSE 
bursts. However, the trend of LGRBs with relatively softer 
spectrum largely weakens in the source frame. 
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